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Abstract 


In this report, we present our findings related to aerodynamic loading of partially inflated balloon 
shapes. This report will consider aerodynamic loading of partially inflated inextensible natural 
shape balloons and some relevant problems in potential flow. For the axisymmetric modeling, we 
modified our Balloon Design Shape Program (BDSP) to handle axisymmetric inextensible ascent 
shapes with aerodynamic loading. For a few simple examples of two dimensional potential flows, 
we used the Matlab PDE Toolbox. In addition, we propose a model for aerodynamic loading 
of strained energy minimizing balloon shapes with lobes. Numerical solutions are presented for 
partially inflated strained balloon shapes with lobes and no aerodynamic loading. 
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Chapter 1 

Aerodynamic loading 


1.1 Aerodynamic loading of natural shape balloons 

1.1.1 Background 

The equations for a natural balloon shape were first derived by researchers at the University of 
Minnesota (see []]). A derivation of these equations is presented in [7, Sec. V]. Extensive com- 
putations were carried out on a digital computer by Smalley (see, e.g., [ 1 5]— [ 1 6]). Axisymmetric 
ascent shapes were considered in [3]. We follow the approach in [3], modeling the balloon as an 
axisymmetric inextensible membrane. However, in the present work, we add caps and aerody- 
namic forces to the model for axisymmetric natural shape balloons. Balloons constructed from a 
single shell, but with sections of different material properties can also be considered. We begin 
with some notation and then provide some physical justification for the governing equations. 

1.1.2 Equilibrium equations 

Since we seek axisymmetric solutions, we need to find a generating curve y — > (z(s),r(s)), where 
y is the arc-length. Let id be the total length of the generating curve. We can express the solution 
as x(s,<t>) (see Figure 1.1), where y € [0,^], <|> 6 [0,2k], and 

*M) = 

ej (<})) = cos<]>i + sin<|>j, 
e 3 ((J)) = k. 

It is useful to define e 2 ((])) to be the unit vector which completes the right-hand triad {ej (4>), e 2 (<]>), e 3 (<(>)}, 

e 2 (<t>) = e 3 (<j>) x ei (<J>) = -sin<t>i + cos<t>j. (1.1) 
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The balloon’s shape is defined by the parametrized surface, 

5={x(5,0) |s€ [0,td],ty€ [0, 2tc] } . (1.2) 

At each point along the curve, s — > x(s, (()), the tangent vector is given by 

3x 

a (5,0) = ^(^,<()) = sin0(5)ei(0)+cos0(5)e 3 (0), (1.3) 

where 0(5) is the angle between a and e3 (t)>) (see Figure 1.1), z!(s) = cos0(s) and t J { s ) = sin0(^). 
The inward pointing normal is 

b(s,4>) = axe2(<l)) = -cos0ej (0) + sin 0e 3 (0). (1.4) 

Let 5o € (0) 4/), and let the vector ni (50, 0) denote the resultant stress defined at each point of 
the curve 0 — > x(so, <}>). Then nj (50, 0) represents the density of the resultant contact force generated 
by the part of 5 with s > so acting on the portion of 5 with s < so- Similarly, let 0 < 0o < 2rc 
and let 112(5, <|>o) denote the resultant stress defined at each point of the curve s — > x(5,0o). The 
vector 112(5, <|>o) represents the density of resultant contact force generated by the part of 5 with 
<t>o < <t> < 0o + £ acting on the part of 5 with <|>o — e < <)) < 0o- We choose e > 0 small enough 
that 0 < <(>o — £ and 4> 0 f £ < 2n. Both ni and n2 are measured per unit length in the deformed 
configuration. 

Because the balloon is modeled as a membrane, we can ignore all bending moments and 
couples. Furthermore, under the assumption of axisymmetry, we can write the contact forces as 

ni(5,<|>) = a m (5)a(5,<J>), 

n 2 (5 , tf>) = o c (5) e 2 (<J>), 

where o m is the meridional stress and o c the circumferential stress. Next, we consider the equilib- 
rium of a test patch A, where 

A = {x(£,y) | 5 0 <£<5, O<\|/<0}, 

for some so e (0,f^) and (|) G (0, 2k). The forces acting on A are the internal forces, ni (5,0) and 
112(5,0), and the external forces, 

f(5,0) = -p(5)b(5,0) - w(5)e 3 (0) -/ D (5)e 3 (0), (1.5) 

where p is hydrostatic pressure due to the difference between the densities of the lifting gas and 
the ambient air, w is balloon film weight, and fo is the aerodynamic pressure density due to drag. 
All external forces are measured per unit area in the deformed configuration. 
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In general, the drag on a submersed body exerted by a fluid that is in motion through the fluid 
is 

F D = ~^pairCpAp\\B ~ V fl/ y|(vg - V air ), (1.6) 

where Cp is the drag coefficient, p a ,y is the density of the air, Ap is the area of the projection of 
the body on the plane perpendicular to the direction of motion, vg is the velocity of the balloon, 
\ a i r is the velocity of the atmosphere (see [7, p. 11-13] and [10, Sec. 4.3]). In our applications, the 
velocity of the balloon will be in the ±k directions. We specify the magnitude of the velocity of 
the balloon |vg| = vg and assume that |v air | = 0. For a descending balloon, Fp is given by 

F d = -\pairCpApv\k. (1.7) 

We define fp{s) — Fp • «3 (<t>), and it follows from Eq. (1.7) that the aerodynamic pressure density 
is 

fD{s) = ~{pairCpV 2 B Ap. (1-8) 

From Eq. (1.8), we see that fp(s) can be expressed as 

fo{s) = - \p a irCpv 2 B \ { _ n/2 fi]{ 9(*)) sinO(s), (1.9) 


where I is the indicator function, i.e., 



x e [a, b] 
xg[a,b] 


We will write fp(Q) to emphasize the dependence of fp on Q = Q(s). The parameter Cp depends 
on a variety of factors (see [7, p. 11-13]); for a descending balloon, we will assume Cp = 0.8. For 
an ascending balloon, fp is similarly defined with Cp = 0.4 and indicator function I(o,n/2]- Note, 
for an ascending balloon fp > 0, while for a descending balloon fp < 0. The caps will be modeled 
as an added thickness, so w depends on s. 


The total force acting on A is 

/'<}> r< l> 

0 = / c m (s)a(s,y)r(s)dy- / a m (j 0 )a(Jo, vM*o)^V 

J 0 Jo 

+ ro.RMd'W- /' o c © 62 ( 0 ) 4 + f S f* t(%, V )r(Qd^d%. (1.10) 

Jsq Jso J $o •'O 

Differentiating Eq. (1.10) with respect to s and (J) and using 9e2/3<)) = — ei, we are lead to the 
following equilibrium equations, 

^(r<7 m a) -G c ei(<|>) + rf= 0. (1.11) 

In this report, when discussing a descending (or ascending balloon), we ignore the dynamics 
of the balloon flight and assume that the balloon is in a quasi-static equilibrium state as given by 
Eq. (1.11). 
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Figure 1.1: Generating curve for a natural-shape balloon (at float) 



1.1.3 The natural-shape model equations 


In this section, we formulate the natural-shape model, including boundary conditions for float and 
descent shapes (ascending shapes can be considered with some slight modifications). Projecting 
Eqs. (1.11) onto the and ei (<{>) directions, we obtain. 


0 = — (ro m cos0) — rw — prsinG — r/b(0), 
as 

0 = — (ro m sin 0) —a c -|- prcosG, 
as 


( 1 . 12 ) 


(M3) 


respectively. After carrying out the differentiation and rearranging terms in (1.1 2)— ( 1 . 1 3), we 
obtain 


0 = -f-(rc m )cos0 — (rc m )sin0^ — rw — prsinG — r/b(0), 
as as 

0 = — (ra m ) sin 0 + (r<j m )cos0— G c + pr cos 9. 

as as 


(1.14) 


If we substitute p — b(z + zo)> where zq is the distance from the zero-pressure level to the bottom 
of the balloon and b is the specific buoyancy of the lifting gas, and solve for dQ/ds and d{r<5 m )/ds, 
we obtain 


(ro m )— = c c cos0- rw sin 0 — br(z + zo)~ f./b(9)sin0, 
ds 

-T-(ro m ) = o c sin0 + rwcos0 + r/o(0)cos0. 


(1.15) 


(1.16) 


While b is assumed to be known as a function of altitude, the value of b depends on a number of 
parameters (e.g., b = b(T a , 7) ; , p a , p ? ) where T a is temperature of the air, T g is the temperature of 
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the gas, p a is the the density of the air, and Pg is the density of the gas). See Appendices A-B for 
further details on the calculation of b. The zero-pressure level za is determined during the solution 
process. Normally , we assume zo = 0 and v# = 0 at float. However, if v# ^ 0, we can also compute 
a natural shape. 

Since the generating curve is parametrized by arc-length, we see that 

~ = cos0 and ^ = sin0. (1.17) 

as as 

It will be convenient to define the total film load T in the meridional direction, i.e., 


T = 2nro m . (1.18) 

Multiplying both Eqs. (1.1 5)-( 1.16) by 2 n, substituting T for 2nra m , and setting o c = 0 (an as- 
sumption of the natural shape), we have 

= — 2nr(ws\nQ + b{z + a) + /b(0) sin0) /T, (1.19) 

= 27tr(w + / D (0))cos0. (1.20) 

For applications considered here, we will assume that a load L is suspended from the base of 
the balloon. Thus, 

7(0) = L/ cos 0 O , (1.21) 

where 0(0) = 0o and 9o is one-half the “cone-angle” at the base of the balloon (see Figure 1.1). 
Note, 9o is not known apriori and must be computed based on certain parameter values and bound- 
ary conditions. If there is no load at the top of the balloon, then 0 = — rc/2 at the top. 


dO 

ds 

dT 

ds 


1.1.4 Problem formulation using the parallel shooting method 

Standard Model 

The governing differential equations for a natural-shape balloon in terms of T are 


dd 

ds 

= — 27tr(wsin0 + 6(z + zo) + r/b(0)sin0)/7\ 

(1.22) 

dT 

ds 

= 27tr(w + /d(0))cos9, 

(1.23) 

dz 

ds 

= COS0, 

(1.24) 

dr 

ds 

= sin0. 

(1.25) 
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At the top of the balloon, we assume 

8(<) = -Jit, 

r(t) = 0. 

The initial conditions for Eqs. (1.22)-(1.25) are 

9(0) = 0 O , 

T(0) = LI cos 0o, 

z(0) = 0 , 

r(0) = 0. 


(1.26) 


(1.27) 


For a given pair (Bo,£) and weight density w(s), Eqs. (1 .22)— (1 .25) can be integrated over the 
interval 0 < s < £ beginning with the initial conditions in Eqs. (1.27). Using a shooting method, 
one determines (©o , £) — (%,d,£d) and the functions 0 (j;0o,^,^), T(s;Bo,d,td)> z(s;&o,d,^d), and 
r(s',Bo,d,£d) that satisfy Eq. (1.22)-(1.25), and Eqs. (1.26), i.e., 


$(&>%, d^d) = -jW, 

r{£d',%,d,£d) = 0. 


(1.28) 


Such a solution automatically satisfies Archimedes principle. For if we multiply Eq. (1.14) by 2tt, 
integrate from s = 0 to s = Id, and used the identity. 




we obtain 


rid 

-T(9 0t d) cos (0 O , d) = 2n wrds -bV-Fo, 

J 0 

and after rearranging terms, 

L + W B = bV + F D , 


(1.29) 

(1.30) 


where Wg is the total weight of the balloon system (shell, caps, and load tapes) and Fp is the total 
drag. 

In the natural-shape model in [3], the weights of the caps and load tapes were added to the 
payload. In the model presented here, the balloon system includes two caps that are modeled as an 
added thickness. We will assume that the number of gores in a complete shape is n g . The weight 
of the load tapes will be incorporated by modifying w appropriately. The lengths of the generators 
for the caps are denoted by ci,C 2 , respectively and are assumed to be known. For a given £, let 
si —£ — C\ and S 2 = £ — C 2 denote the location of caps’ edges. This means that w(s, -(-0) ^ w(s/ — 0) 
for i — 1,2 and dB/ds and dT jds are discontinuous (another discontinuity arises when fo 0). 

Next, we outline a parallel shooting method for computing a design shape with a variable 
thickness. See [13] and [14] for a discussion of the parallel shooting method. In the following, we 
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denote the length of the first, second, and third interval by l\ = s\, I 2 = si- s\, and It, = t — S 2 , 
respectively. We divide n g w, by 2nr(s) in order to get the average weight density of tape with 
respect to area. The total length t is not known beforehand and will be computed as part of the 
solution process. We define w(s) as follows, 

w(s) = Wfji m {s)+n g w t /(2nr(s)), 0 <s<£, 


where 

{ Wf, 0 < s < l\ 

2 -Wf, l) <s<l 2 , (1.31) 

3w/, I 2 < s < £ 

Wf is the weight density of a single layer of balloon film, w t is the load tape weight density and n g 
is the number of gores. Eq. (1.31) can be easily modified to accommodate balloons with varying 
material properties, other than the thickness due to the caps. Eqs. ( 1 .22)— ( 1 .25) are equivalent to 


dQ 

—2 = -2 Tin ( Wi sin 0/ + b{zi + zo) + foi®;) sin 0;) /7}, 
ds 

— = 2%ri (w, + /o(9/) ) cos 0 / , 
ds 

dzi 
ds 
dri 
ds 


- cos 0 ,-, 

= sin 0 ,, 


(1.32) 

(1.33) 

(1.34) 

(1.35) 


for Si- 1 < s < si and i= 1,2,3, with the convention that .so = 0 and 53 = t. For a function such 
as w(s), we will follow the convention that Wi(s ) = w(s) for j,_i < s < si, and similar defini- 
tions for the other relevant functions 0, , Ti,n, r,-. Eqs. (1.32)-(1 .35) are supplemented by the initial 
conditions. 


01 (0) = 

«?. 


T\ (0) = 

L/ cosG®, 


n (0) = 

0, 


21 (0) = 

0, 


0i(si) = 

8?+|. 

1 = 1)2, 

Tifa) = 

7"0 

i= 1,2, 

Zi(S') = 

7° 

z i+ 1’ 

f = 1,2, 

n(si) = 

r° 

r i+ 1’ 

i= 1,2. 


(1.36) 
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Eqs. (1.36) involve nine explicit free parameters {©pG®, 7j°, T^z^z^f r 2 } anc * * m P^ clt P a_ 
rameter i. Using the conditions at the top of the balloon given by Eq. (1.28), we are led to ten 
auxiliary conditions involving the unknown parameters, 


Oi(si) = 

9?+|. 

1=1,2, 

Ti(sd = 

l£|. 

1= 1,2, 

n(st) = 


* = 1,2, 

Zi(si) = 

7 ° 

z i+ 1 » 

*=1,2, 

r*(0 = 

0 , 


e 3 W = 

- 5 *- 



By construction, the functions 0(s) = 0,-(s), T (s) = 7}(s), r(s ) = n(s), and z(s) = Zi(s) for * = 1,2,3 
are continuous on 0 <s <ld- 


Following the approach of [13], a new independent variable s is introduced on each interval 
by the transformation 


s= (s-Ti)/li, i= 1,2,3, 


where Ti = 0,t2 = l\ and T3 = l\ + h- By construction, 0 < s < 1 on each subinterval. Let z,i = 
n{s ) , n = n(s) . . ., and so on. In the rescaled coordinates, we find Eqs. ( 1 .32)-( 1 .35) can be written 
as 


dfy 

ds 

dTi 

ds 

dzi 

~d§ 

dr\ 

ds 


= —InfiU (w,- sin0/ + b(zi+zo) + /£>(0/) sin0,) /7}, 
= 2nnli (wi + /b(0i)) cosG/, 

= /,COS0 M 

= // sin 0 ( -, 


(1.38) 

(1.39) 

(1.40) 

(1.41) 


for t = 1,2,3. 
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The auxiliary conditions in the rescaled coordinates are 


0/(1) = 

8?+l> 

i = 1)2, 

U 1) = 

T'O 

i i+l> 

f = 1,2, 

h(\) = 

r° 

r i+ V 

1=1,2, 

*(i) = 

7 ° 

z i+ 1’ 

i= 1,2, 

M 1 ) = 

o, 


03 (*) = 



2 n. 



The problem of computing a standard natural-shape balloon is to find a solution of Eqs. (1.38)- 
(1.41) on the interval 0 < f < 1 satisfying Eqs. (1.42) with vg = 0 and zo = 0. 

Natural shape design with drag 

Eqn. (1.30) suggests there is a one-parameter family of design shapes, parametrized by vg. 
We make the standard assumptions that c c = 0 and zo = 0, but do not assume that vg is zero. The 
net effect is that the volume of the design shape is reduced and the balloon is taller and narrower 
for decreasing vg. Natural shape balloons with drag will be computed in Section 1.5.2. A natural 
shape design with drag is a solution of Eqs. (1.38)-(1.41) on the interval 0 < s < 1 satisfying 
Eqs. ( 1 .42) with vg ^ 0 and zo = 0. 


1.1.5 Partially inflated natural shape balloons 

In the following, we assume that the design shape is known and the generating curve for the design 
shape is 

(R d (s),Z d (s)), 0<s<e d . 

If vg = 0, in our exposition, we assume that 

bV = b d V d (1.43) 

where V d is the design volume, b d is the buoyancy at float, and V is the current volume. Eq. (1.43) 
implies that b = b d V d / V. Since V < V d , the balloon envelope is not fully deployed. If vg ^ 0, 
we assume b = xbk, and Archimedes Principle and Eq. (1.30) imply that the volume must satisfy 
V=(W B +L-F D )/b. 
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To account for excess material, we modify the film weight density (replacing w(.y) — > w(s)Rd(s) /r{s)). 
In particular, Eqs. (1.22)-(1.25) are replaced by (see [3]) 


d9 

ds 
a IT 

ds 

dz 

ds 

dr 

ds 


- 2n(Rd (5) w sin 0 + br(z + zo) + rfD (0) r sin 0) / T, 
2n (R d (5) w + rf D (6) ) cos 0, 

COS0, 

sin0. 


( 1 . 44 ) 

( 1 . 45 ) 

( 1 . 46 ) 

( 1 . 47 ) 


Since we assume that the generating curve is inextensible, id is constant, and we introduce 
the zero-pressure level zo as the parameter to be determined. In particular, we replace Eqs. (1.26), 
by 

= -3*, 


(1.48) 


0) = 0. 


When the volume of the balloon gets very small, 0® is nearly zero, and it virtually impossible 
to solve Eqs. (1.48) with an ordinary shooting method. This problem was circumvented in [3], by 
subdividing 0 < s < id and solving the governing equations with a parallel shooting method. For 
the present application, we subdivide the three intervals [0,/t], [h,h], [h,id\- For a given subdivi- 
sion obtained by adding adding new break point Sk, four additional parameters {0°, 7^°, r^z®} are 
introduced with four additional auxiliary equations, 

eS-ifa) = ^ 


( 1 . 49 ) 

rlM) = % 


Z k - 1 — 


7 ° 

Z k- 


After relabeling of quantities, for n subdivisions there will be a total of 2 + 4(n — 1 ) parameters to 
be determined, i.e., 


_ aO aO 7 -O 0 _0 aO 

20? y l I 2 i r 2i z 2i • • 


T’O 0 0 

*n—\ ’ r n - 1 ’ z n — 1 • 


The auxiliary equations are rescaled in the same manner as Eqs. (1.38)-(1.41). 
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1.1.6 Numerical results 


Standard design shape 

In the following, we present results on descending axisymmetric inextensible natural shape bal- 
loons using the parallel shooting method described in Section 1.4. We consider a design shape 
based on the parameters in Table 1.1. The values in Table 1.1 correspond roughly to a solar Mont- 
golfiere balloon made of 3.5 micron mylar with a payload of 6 kg. The examples were chosen 
to illustrate how the shapes varies as certain parameters change. We first compute a design shape 
with vb = 0. The height and diameter of the design shape is shown in Table 1.1. We then compute 
shapes for x = 2 and 1 = 5. When vg = 0 and x = 2, we find V = 0.5 V<j. When vg = 0 and x = 5, 
we find V = 0.2 Vj. In Figures 1 .2- 1 .3, the design shape profile is shown to the left of the partially 
inflated ascent shape. Plots of 0 (j), r(.y),r(.s),x(s) are also shown. 

In Figure 1.2(a), we present an ascent shape with V = 0.5 Vj. To eliminate drag from our 
model, we set Co = 0. The variation in shading indicate the subintervals used in the parallel 
shooting method (8 subintervals were used). Next, we set Co = 0.8 and computed a corresponding 
solution shown in Figure 1.2(b). Because of Archimedes Principal (see Eq. (1.30)), if we maintain 
the same b and increase F D , then the effective volume V must decrease. In the zero drag case, 
the total meridional tension increases monotonically. With drag acting on the lower portion of the 
balloon, the tension is reduced. Next, we increased the velocity and air density and considered 
shapes with a smaller volume. In Figure 1 .3(a), we present an ascent shape with V = 0.2 Vj. Next, 
we set Co = 0.8 and computed a corresponding solution (see Figure 1.3(b)). In Table 1.2, we 
compare quantities related to the shapes presented in Figure 1 .2 and Figure 1.3. 

Remark. It should be noted that the assumption, b V = bj Vj is more relevant to ascent 
shapes than for a descending balloon that is being inflated. For a descending balloon that is being 
inflated (such as a Solar Montgolfiere), the higher buoyancy values occur when the balloon is near 
full inflation and the low values occur for low volumes. Nevertheless, the results presented in this 
section are for demonstration purposes and illustrate how such shapes can be analyzed. 

Natural shape with drag 

In this section, we consider a balloon descending with a constant velocity and compute the cor- 
responding design shape. We assume a balloon with a 6 kg payload and a balloon envelope with 
density 8 kg/m 2 . Additional design parameters are presented in Table 1.3. 

In Figure 1.4, we present a family of design shapes parametrized by vg. In Figure 1.4(a), 
vb = 0. In Figures 1.4(b)-(d), we present design shapes with velocities -2 m/sec, -4 m/sec, and -6 
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Table 1.1: Design Parameters for Test Case 


Parameter 

Description 

Design Value 

bd 

Buoyancy 

3.9020375e-04 lbf/ft 3 

L 

Payload 

13.22750 Jbf 

Wf 

Film weight density 

1.43369e-03 lbf/rf 

W, 

Tape weight density 

0 lbf/ft 

n 8 

Number of gores 

27 

v d 

Volume 

62290 ft 3 

e d 

Gore length 

76.738 ft 

0o 

Angle at base 

59.32 deg 

maxZj 

Balloon Height 

45.76 ft 

2 • max Rd 

Diameter 

53.219 ft 

At 

Surface Area 

7727.26 ft 2 

W T 

Gross Weight 

24.306 lbf 


Table 1.2: Comparison of descending shapes with and without drag 


Quantity 

Description 

0.5 Vd 

0.2 V d 

C D 

Drag coefficient 

0 

0.8 

0 

0.8 

vq (ft/sec) 

Velocity of balloon 

0 

-4.59 

0 

-16.40 

F d (lbf) 

Total drag 

0 

7.06 

0 

21.91 

zo (ft) 

Zero-pressure level 

29.97 

34.47 

42.48 

57.37 

b (lbf/ft 3 ) 

Buoyancy 

2 b d 

2 b d 

5b d 

5b d 

max /-(ft) 

Max radius 

19.86 

17.60 

14.40 

6.64 

maxz (ft) 

Height 

58.11 

60.74 

63.69 

71.99 
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0.0007804 tau: 2 tau = 0.500 z = -34.7 F n : 7.0591b ncalls: 151 


(a) Without drag 


C D : 0.0 Rho -r = 0.0611 kg/m 3 v: -1.40 m/s 0vsS 



(b) With drag 



Figure 1 .2: Comparison of descent shapes with and without drag 











0.001951 tau: 5 tau = 0.200 F • -57.371 lb ncalls: 2.191092e+01 b = 227 tau 

















Table 1 .3: Parameter values for a natural shape design with drag 


Quantity 

Value 

Payload 

6 kg 

Film density 

8 g/m 2 

Pair 

0.00827 kg/m 3 

P gas 

0.00541 kg/m 3 

b 

280.764 kN/m 3 


Table 1.4: Comparison of natural shapes with drag as a function of velocity 


Velocity (m/sec) 

0 

-2 

-4 

-10 

Gore length (m) 

36.87 

30.08 

33.86 

23.68 

Go (deg.) 

68.06 

67.03 

63.84 

43.60 

Height (m) 

20.42 

20.06 

19.07 

14.88 

Diameter (m) 

26.63 

26.05 

24.38 

16.12 

Volume (m 3 ) 

7291.98 

6811.29 

5555.57 

1619.63 

Surface Area (m 2 ) 

1856.87 

1774.93 

1551.50 

689.73 

Mass (kg, Mb + M^ d) 

20.85 

20.20 

18.41 

11.52 


m/sec. As the velocity increases in magnitude, the volume of the design shape is reduced and the 
ratio of diameter to height is also reduced. Once a design shape has been computed, we can also 
determine partially inflated natural shapes as we did in the previous section. 

In Table 1.4. we present computed quantities related to the family of design shapes. 

1.1.7 Remarks on natural-shape analysis 

It is interesting to compare the results from Table 1.4 with that of a sphere of comparable volume. 
For example, consider the case when vb = —10 m/sec. In this case, the volume, diameter, and 
balloon weight of the natural shape is 1619 m 3 , 16.12 m, and 5.517 kg. For a sphere of the same 
volume, the diameter is 14.56 m and its weight is 5.334 kg. Even thought the natural shape weighs 
only 0.18 kg more than the sphere, it produces a drag equivalent to 6.89 kg, versus 5.623 kg for 
the sphere. The natural shape produces over 22% more drag than a sphere with the same volume. 

Aerodynamic drag can be incorporated into a model for inextensible partially inflated natural 
shape balloons. While real balloons are not inextensible and the hoop stresses are not zero, the 
results in this section can be applied to elastic balloons that are not axisysmmetric. In Section 2, 
we show how to extend this approach to our EM-model for a more general class of balloon shapes. 
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1.2 Aerodynamic loading of partially inflated balloons with lobes 


In this section, we outline a mathematical model for the aerodynamic loading of a lobed balloon 
shape. For simplicity and ease of exposition, we focus primarily on the case of drag on a descend- 
ing balloon. We adapt the method of Chapter 1 to balloon shapes with and without lobes. Our 
approach is to derive an expression for the potential energy due to drag, so it can be added to our 
model for lobed energy minimizing shapes (i.e., LEM-shapes). 

1.2.1 Background 

A feature characterizing partially inflated scientific balloons is a lobe structure surrounding the 
gas bubble. When the gas bubble is of insufficient volume, the balloon envelope will collapse, 
forming folds of excess material and lobes. We present a mathematical model for a balloon shape 
with such features that is subjected to aerodynamic loading. We formulate a variational principal 
for the balloon system, and seek to determine a shape that minimizes energy subject to a volume 
constraint. The energy includes film and load tape strain, film and load tape weight, hydrostatic 
pressure, and drag. 

Based on observations of real balloons, it is reasonable to assume that lobed structures are 
periodic. Balloon shapes with lobes were investigated in [5] using variational techniques where 
the energy of the balloon system included hydrostatic pressure and film weight. The energy of 
the balloon system was minimized subject to constraints on the volume and the lengths of certain 
fibers. Geometrically, the shapes in [5] possessed many features seen in real balloons: lobes 
surrounding the gas bubble, a sphere-like top, internal folds of excess balloon fabric, and flat wing- 
like structures (fins) hanging beneath the gas bubble. However, the formulation of the problem did 
not include strain energy. In [2] and [4], the film and load tape strain energies were considered for 
shapes modeled by a half-gore. In this work, we formulate a variational principal for a partially 
inflated balloon subject to aerodynamic drag in such a way that a lobed structure can (but not 
necessarily) form. While drag influences the balloon shape, lobes will form independently of drag. 

In the following, let S(V,q,r) denotes a deformed balloon with volume V, q lobes, and r 
gores per lobe. The total number of gores is n g and n g = rq. If it is clear from context, we will 
omit certain arguments and write 5 instead of S( V,q,r). Sp denotes a fundamental section of 5 
(Sp is centered about a vertical plane of reflectional symmetry). By symmetry arguments, r is odd, 
i.e., r = 2k + 1 . A superscript j = 1 , 2 , . . . , r denotes a quantity related to the y'th gore of Sp. When 
r = 1 , the balloon shape can be modeled by a half-gore. When r > 1 , the balloon can have lobes. 
Viewed from above, the discrete design shape will have the symmetries of a regular polygon with 
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n g sides. In general, 5( V,q,r) will be invariant under the dihedral group D q and viewed from 
above, it will have the symmetries of a regular polygon with q sides. To simplify our exposition, 
we will consider a natural-shape design with n g = 27 (note. On C Cg C 5?). A shape with Dg 
symmetry has nine lobes and each complete lobe is constructed from 3 gores. 


1.2.2 Potential energy for aerodynamic drag 


Aerodynamic drag can be introduced into our balloon shape model by adapting Eq. (1.6) appro- 
priately. Consider an arbitrary facet in the triangulation of a balloon shape that is not part of a 
fold. Let vo, vi , V 2 denote the vertices, and eo = v\— vo, ?i = V 2 — V[ , ?2 = vo — V 2 . We assume that 
v 0 , vi, V 2 are ordered in such a way that the unit normal n = <?o x e\/\eQ x e \ | points outward. For 
a balloon descending with a velocity v B = — Vflk, the drag on typical triangle T is: 



2 Piur('l) v jj dA • k, 
0 


0 < arccos(— k • n) < 7t/2, 
otherwise, 


(1.50) 


where dA = ^\eo x e\ |n. Adding up the drag contributions of each facet that contacts the atmo- 
sphere yields the total potential energy of the drag on the balloon, i.e., 

E D {S) = -\f s PairC D v 2 B dA'k. (1.51) 


1.2.3 Mathematical model for aerodynamically loaded balloon shapes 


In Table 1 .5, we present design parameters that are used for the illustrations in this section. In the 
following, the reader should refer to Figure 1 .5 for notation related to the reference configuration. 
Figures 1 .6(a)-(b) shows a 3 gore lobe of a 27 gore design and Figures 1 .7 shows the complete 
shape. 


Let V; , V mt , V m : , V m , V mr , V r denote vertices along a circumferential fiber at the ith station of 
Gore-y in the flat reference configuration. Let v/,v m /,v m «,v m ,v mr , v r be the corresonding vertices 
in the deformed configuration. Note that V m i and V m are the same material point. Vertices v m >,v m 
represent the end location of the respective left-half and right-half circumferential fibers (if these 
fibers were to be extended in a straight line). Vertices of the type v mr and v r have three degrees of 
freedom (DOF). Vertices of the type v m i and v/ have 0 DOF, since v m / is equal to v mr on Gore-y 
and v/ in Gore-y is equivalent to v r in Gore-y'-f 1. Vertex V mr — (Y mr ,Z mr ) has one DOF, since Y mr 
is allowed to vary between Y m and Y r . Similarly, V m i has 1 DOF. Vertices in the form v m < and v m 
are determined as follows: 

V mri Ymr ^ 

v m — Vmr (ymr v r) ' ^(V/n, V mr , V r ), Y m <C.Y mr <CY r , 

(-X/nr (Xr COs(7l/ q),ymr (Y r Y m ) sin(j l/q),Zmr) , Y mr ^ Y r , 
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Table 1.5: Design Parameters for a 27 gore balloon 


Parameter 

Description 

Design Value 

b d 

Buoyancy 

1.9459e-02 lbf/ft 3 

L 

Payload 

107 lbf 

Wf 

Film weight density 

0.00384 lbf/ft 2 

Wt 

Tape weight density 

0.00301 lbf/ft 

n g 

Number of gores 

27 

Vd 

Volume 

5968.46 ft 3 

td 

Gore length 

36.0028 ft 

e 0 

Angle at base 

51.44 deg 

maxZj 

Balloon Height 

22.935 ft 

2 • maxRj 

Diameter 

23.764 ft 

At 

Surface Area 

1617.29 ft 2 

W t 

Gross Weight 

116.136 lbf 


Figure 1.5: A segment of a gore in the reference configuration. 
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Figure 1.6: LEM-shape at V = 0.15 Vj for a 27 gore design, (a) Inside view of 3-gore lobe; (b) Outside 
view of a 3-gore lobe. 


Figure 1.7: A complete LEM-shape at V = 0.15 V d for a 27 gore design with nine lobes. 
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v m i 


{xml (Yt ^m') cos(lt/q),y m t (Y[ Y m > ) sin(jc/q) , Zml) i Y m > Y/, 

' + (v m ; V/) • X(V m ', V m /, V/), T/ < Y m i < Y m >, 


where T(V m ,V mr ,V r ) — (Y mr — Y m ) / (Y r — Y mr ), and x(V r m /, , V/) — (1^,/ Y m i)/(Y r Y mr ). The 
vertex that lies at the top of the balloon has one DOF, since it can vary along the z axis only. Linear 
equality constraints for the bounding planes of a fundamental section Sf are in the form: 


ymr = o, (1.52) 

y r — ian{K/q)x r = 0. (1.53) 

For lobed shapes, we impose the condition, 0 < y < tan(7t /q)x, for vertices (x,y,z) with non-zero 
degrees of freedom that are not constrained by Eqs. (1.52)-(1.53). These conditions insure that a 
balloon shape will have /^-symmetry. The variational principal for an energy minimizing balloon 
shape is to minimize the total energy of the balloon system subject to a volume constraint over a 
class of shapes C q that incorporates geometry, symmetry, folds, etc. Next, we write out explicitly 
the terms in the variational principal. To simplify notation, we use S to denote both the actual 
surface and its approximation. The /th triangle of the yth gore in the reference configuration is 
denoted by T t J and the corresponding triangle in the deformed configuration is denoted by T t J . The 
reference configuration is Sf = U^ =1 U^, T( and the deformed configuration is Sf = U^_, U^j T t J . 
The total number of facets in a typical gore is N. The total number of facets in a typical gore that 
lie on the outside of the balloon is N°. The current target volume is V. Components of the total 
energy of the balloon system contributed by the yth gore of Sf are: 


(a) Hydrostatic Pressure: 



(b) Film Weight: 

4ta = */£^“W)area(r/) t 

/= 1 

(c) Tape Weight: 

n c + 1 . 

E { tapes w t ^ Zj,2 \E[ I’ 

i=l 

(d) Top Fitting Weight: 

Etop — WtopZtopi 


(g) Aerodynamic Drag: 

E, drag ~ ~ 2 9airCoV[) ^ area ^Pxy[^;/] j , 
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where P X y[ T] is the projection of facet T onto the xy-plane and the summation in l' is over all 
facets with unit normal making an angle less than n/2 with the — k direction (see Eq. (A.3), 


(e) Tape Strain: 


(f) Film Strain: 


1 

Stapes = £ (e/) 2 |£/|, 


1= l 


AT 


s{ita = 2 2 (0(7-/)n(j;') : r(r/)area(3-/). 
/=1 


The volume is W = 2 V/ 7 , where n c is the number of circumferential fibers, wy is the film 

weight density, co(T) is the number of film layers in facet T, zj is the centroid of T/, w, is the tape 
weight density, K, is the tape stiffness, n is the Second Piola-Kirchoff stress, y is the Green strain, 
e/ is the strain in the *th segment of the yth load tape, z/ i 2 is the centroid of the ith meridional 
segment in the yth load tape, | Ej | is the length of the ith segment of the yth load tape, b is the 
buoyancy, V( is the volume contribution due to facet l of Gore- j. The height of the top fitting 
above the base of the balloon is Zt 0 p and its weight is w lop . Summing the contributions over all 
gores within a lobe and multiplying by q, we obtain the energy of the complete balloon: 


r_ 


(1.54) 


and the volume constraint: V = <?Sy=i V*- When r — 1, we call a shape S that minimizes Ej ota [ 
an energy minimizing shape (EM-shape). When r > 1, the shape has lobes and we refer to it as a 
lobed EM-shape (LEM-shape). 


1.3 Two-dimensional potential flows 

In this section, we consider two dimensional uniform flow around an axisymmetric balloon shape. 
Bernoulli’s Law can be used to compute a pressure distribution. 

1.3.1 Potential flow and the PDE Toolbox 

In Fig. 1.8(a), we present the profile of a natural shape balloon (near full inflation) and in Fig. 1.8(b), 
we present a partially inflated natural shape (at roughly 5%-inflation). We assume that the balloon 
is stationary and far away from the balloon, the fluid moves parallel to the x-axis at a speed of 
about 4 m/sec 
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(a) 


(b) 


Figure 1.8: Axisymmetric balloon profiles: (a) full inflation; (b) 5% inflation 


Let B be the rectangle, in the xy-plane, 

B = {(*,y) | L\ < y < L 3 , L 2 < x < L 4 } 
and r)B = dB\ U dB 2 U dB 3 U dB 4 , where 

= {(x,y) GB,y = Li}, 
dB 2 = {{x,y) EB,x = L 2 }, 
dB 3 - {{x,y)eB,y = L 3 }, 

= {(jc,j) G B,x = L 4 }, 

where L\ + L 3 = 0. 

Let A denote the interior of the balloon shape (see Fig. 1.8), and let Q be the region exterior 
to A, but interior to B: 

Q = B\A. 

In Fig. 1.9(a) and Fig. 1.9(c), we present a discretization of the region Q produced using 
the PDE Toolbox that is exterior to the balloon for full and 5% inflation, respectively. The initial 
meshes can be refined by the user or adaptively refined by PDE Toolbox (see Fig. 1.9(b) and (d)). 
For the purpose of this exposition, we impose Neumann boundary conditions on dB 2 UdB 4 , and 
Dirichlet boundary conditions on 3A and dB\ U dB 3 . We are led to the following boundary value 
problem: 

Av|t = 0, on Q, 

\|/ = 0, on 3A, 

(BVP) = 0, ond/?2Udfi 4 , 

an 
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Figure 1.9: From PDE Toolbox (a) Initial mesh for full inflation shape; (b) Adapted/refined mesh 
for full inflation shape; (c) Initial mesh for 5% inflation shape; (d) Adapted/refined mesh for 5 % 
inflation shape. 
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VI \f = -V|/o, on dBi , 
i|/ = \j/ 0 , ondBj, 


where \|/ is the stream function, and d/'dn is the normal derivative. In this case, the velocity field is 


v = 


a\|/, a\|/. 


(1.56) 


The PDE Toolbox can be used to solve BVP, compute the velocity field, and to plot the corre- 
sponding flow lines on the adapted mesh. The results are shown in Fig. 1.10. Under the assumption 
of steady irrotational flow of a fluid with uniform density p, we have the following 


- + f |v| 2 = Constant, 


(1.57) 


from which we can determine the pressure P\ v is the velocity field (see Eq. 1.56). Eqn. (1.57) is 
called Bernoulli’s Law. 


1.3.2 Concluding remarks 

The PDE Toolbox is restricted to two dimensional flows, and we would have difficulties extending 
this approach to a flow problem for a balloon with lobes. As outlined in Chapter 2, the methods in 
Chapter 1 can be extended to full three dimensional problems. Drag is most likely the dominant 
aerodynamic load for the inflation/deploymentphase of a planetary balloon mission and so it makes 
sense to develop further the methods of Chapters 1-2. 
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1.4 Dynamic model for a descending spherical balloon 


In the following, we consider a simple dynamic model for a descending spherical balloon that is 
inflating with air. The motivation is a Mars solar Montgolfiere balloon. The total force on the 
balloon is equal to the sum of the lift due to the gas, dead weight of the balloon system (payload 
and balloon), and drag. In particular, we have (see [7, p. 11-15, Eq.(32)]) 


d 2 z 


Fuft ~g(mB+m L )+F D , 


0<t<t Rd 


(1.58) 


where 


z 

t 

m B 

m L 

r 0 

m v 

V 

Pa 

R 

Ad 

T a 

m g 

P g 

Fuft 


Fd 

Cd 


z(t) (altitude of balloon, geopotential height) 
time 

mass of balloon 
mass of payload 

radius of intake at base of balloon 

m B + ntL + m g + C m pa(z) V (virtual mass) 

Kr o (z(0 - zo) + Vq(Vq is the initial fill) 
p a (z) (density of atmosphere) 

3K\ 1/3 


( 


47t ) 


(radius of sphere with volume V) 


= jimax(/? 2 , r 2 ) (frontal area of spherical balloon) 

T a (z) (temperature of atmosphere) 

7 f Z ° 

m gt o + nr~ J p a {H)dH 

m g f V (mean density of gas) 

gp a V-gm g 

gV p a 


Q 


r a +e 


Jp a{z)Co 


dz 

Jt 


dz 

dt 


Ad 


0.8 (for a descending balloon) 

time when the balloon first reaches full inflation, i.e., R = Rj 
Units: Length-meters, Mass-kilograms, Time-seconds 
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Table 1 .6: Parameter values for a spherical balloon 


Quantity 

parameter 

Value 

Mass of payload 

m L 

6 kg 

Diameter (at full inflation) 

2 Rd 

15 m 

Mass of balloon 


5.65 kg 

Film density 

P / 

8 g/m 2 

Initial velocity 

VB, 0 

-50 m/sec 

Supertemperature 

O 

112° K 


Table 1 .7: Test Case data 


Quantity 

parameter 

Value 


t 

z 

z' 

z" 

0 

36520 

-50 

-9.230 

87.6 

34957 

-10.99 

0.068 

720 

30501 

-4.06 

0.008 

1200 

29378 

-0.712 

0.007 

2400 

29356 

-0.030 

-0.001 


The term C m p a (z) V is the added environmental mass that is assumed to be accelerated along 
with the balloon, payload, and gas masses. For balloon applications, one typically assumes C m = j 
([7, p.g 11-12]). The quantities p a (see Eq. (C.5)) and T a (seeEq. (C.l)) are discussed in Appendices 
A-B. For t > tR d , we assume that the balloon maintains a spherical shape with a radius of R = Rj, 
and terms in Eq. (1 .58) can be appropriately modified. 

Next, we consider a test case based on the values presented in Table 1.6. When released, 
we assume that the balloon is falling with a velocity of —50 m/sec. At the altitude where the 
balloon is fuly inflated, we assume that the atmospheric pressure is 0.00827 mb; this corresponds 
to an altitude of 34957 m. We estimate that the balloon must fall 1563 m in order to reach full 
inflation, and so we assume the release altitude is 36520 m. Integrating the Eqn. (1.58) with the 
iniital conditions of z'(0) = —50 and z(0) = 36520, we find that it takes 87.6 sec for the falling 
balloon to inflate, i.e. tR d = 87.6). After full inflation, ^(87.6) = —10.99 m/sec., z"(87.6) = 0.068 
m/sec 2 . Other values of z,zf,z" are displayed in Table 1 .7. After about 720 sec, the balloon reaches 
is minimum altitude and begins slowly oscillating about some mean altitude. In Figure 1.11, we 
present data related to the U.S. Standard Atmosphere. In Figure 1.12, we present plots of (t,z(t)) 
and {t,z'(t)) for 0 < t < tR D . 
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Figure 1.12: Altitude and velocity from test case; (a) z{t) - altitude; (b) { ( t ) - velocity. 

1.5 General Remarks on Aerodynamic Loading 

The methods outlined in Chapters 1-2 are feasible approaches to incorporating aerodynamic load- 
ing into balloon shape analysis. While the natural shape model described in Chapter 1 does not 
include elasticity, it provides a good first step for determining the effects of aerodynamic loading. 
The approach outlined in Chapter 2 can be adapted to nonaxi symmetric shapes and fits very nat- 
urally in our program for modeling partially inflated balloon shapes by variational techniques and 
energy minimization. 

Adding these aerodynamic capabilities to our LEM-model, we would be able to provide an- 
alytical capabilities that could be used to support missions such as the Mars Solar Montgolfiere 
Balloon that is being studied by JPL, especially for the critical inflation/deployment phase. These 
capabilities could also be applied to ultra long duration flights, including in-spool configurations, 
constrained ascent shapes, and ascent shapes. While the analytical results of Chapter 1 are related 
to axisymmetric shapes, with capabilities that are now being developed in LEMsolver and the ap- 
proach outlined in Chapter 2, we should be able to develop a tool for analyzing aerodynamically 
loaded balloons for a variety of configurations. 
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Chapter 2 

Lobed Balloon Shapes 


A feature characterizing partially inflated scientific balloons is a lobe structure surrounding the gas 
bubble. We present a mathematical model for such a balloon shape. We formulate a variational 
principal for the balloon system, and seek to determine a shape that minimizes energy subject to a 
volume constraint. The energy includes film and load tape strain, film and load tape weight, and 
hydrostatic pressure. Numerical solutions are presented for a 24-gore balloon with and without 
lobes. 


2.1 Introduction 

Ascent shapes of large scientific balloons are characterized by the formation of lobes surrounding 
the gas bubble. Based on observations of real balloons, it is reasonable to assume that these struc- 
tures are periodic. Balloon shapes with lobes were investigated in [5] using variational techniques 
where the energy of the balloon system included hydrostatic pressure and film weight. The energy 
of the balloon system was minimized subject to constraints on the volume and the lengths of cer- 
tain fibers. Geometrically, the shapes in [5] possessed many features seen in real balloons: lobes 
surrounding the gas bubble, a sphere-like top, internal folds of excess balloon fabric, and flat wing- 
like structures (fins) hanging beneath the gas bubble. However, the formulation of the problem did 
not include strain energy. In [2] and [4], the film and load tape strain energies were considered 
for shapes modeled by a half-gore. In the present work, we formulate a variational principal for a 
partially inflated balloon in such a way that a lobed structure can (but not necessarily) form. Based 
on observations, it appears that a lobed shape is somehow preferred over a cyclic shape without 
lobes when the volume is significantly less than the float volume. 

Applying our model to a representative shape (a 24-gore balloon design), we investigate by 
parametric studies the existence of a minimum energy selection principal behind the formation of 
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lobes in balloon shapes. In addition, we compare various quantities such as total energy, strain 
energy, and maximum principal stress resultants for balloon shapes with a different number of 
lobes. When there are a large number of gores per lobe, we found that shapes with self-intersections 
could evolve. By monitoring the evolution of a family of LEM-shapes and (when necessary) adding 
a perturbation to break the symmetry, we found that in the cases studied, the self-intersection could 
be avoided. 


In the following, we let S( V,q) denotes a deformed balloon with volume V and q lobes. The 
number of gores per lobe is p and the total number of gores is n g = pq. If it is clear from context, 
we will omit certain arguments and write S instead of S( V,q). Sp denotes a fundamental section 
of S ( Sf is centered about a vertical plane of reflectional symmetry). In the past, we assumed 
that p was odd. However, in this treatment, we will only require that q divides n g . A superscript 
j— 1 , 2 , . . . , p will be used to denote a quantity that is related to the y'th gore of Sf- When p = 1 , 
the balloon shape can be modeled by a half-gore. When p > 1 , the balloon can have lobes. Viewed 
from above, a design shape will have the symmetries of a regular polygon with n g sides. In general, 
S( V,q) will be invariant under the dihedral group D q . Viewed from above, S( V,q ) will have the 
symmetries of a regular polygon with q sides. The set of balloon shapes invariant under D q will 
be denoted by C q . To demonstrate our approach, we will consider a natural-shape design with 
n g = 24, noting 1 

C 2 4 C 2 C C 6 C C-i and C 24 C Cg CC C 4 . (2.1) 

We compute a family of balloon shapes with D 24 symmetry and a family of balloon shapes with 
Dk symmetry and k = 24, 12, 8, 6, 4, 3 The numerical solutions we present here need not be ail the 
possible solutions. In Table 2.1, we list the combinations of p and q that are possible for a 24-gore 
balloon. Table 2.1 contains design shape parameter values. Before we present our mathematical 
model, it is a good to have a rough idea of the geometry of a partially inflated shape for various 
inflation levels. Beginning at float altitude, we assume V = V d (its design volume). At float, 
the balloon envelope is fully deployed. As the volume decreases, the balloon becomes taller and 
narrower. Shapes of this kind can be modeled by a half-gore and were the focus of [4], For 
0.6 K, < V < Vd, we see that as the volume decreases, each deformed gore moves toward the 
z-axis and an internal fold of excess material forms in the center of the gore. In an internal fold, 
the exterior of the balloon comes into contact with itself. To reduce the number of unknowns, in 
our treatment here, we will not use an explicit model of the internal fold. For a sufficiently small 
volume, the inside of the balloon will come into contact with itself. Since the precise nature of the 

1 In the variational formulation of the balloon problem, the energy of the balloon shape will be minimized over some 
C q . However, due to Eq. (2.1)), it is possible in certain cases that a shape 5^/ with > q to evolve when minimizing 
over Sq. If this does happen, we can add a perturbation to an intermediate shape in order to break the D symmetry. 
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Gores per Lobe 

P 

Number of Lobes 
<7 

1 

24 

2 

12 

3 

8 

4 

6 

6 

4 

8 

3 

12 

2 


Table 2. 1 : Parameter values for balloon 


E (psi) v 

Wf (lbf/in 2 ) 

w, (lbf/in) 

K, (lbf) 

e (in) 

b d (lbf/in 3 ) 

K/(in J ) 

n g £ d (in) w wp (IbO 

36000 0.82 

2.667e-05 

2.50833e-04 

2810 

0.0008 

1.1330556e-05 

10775165.571 

24 433.15 1.5 


distribution of folded material is unknown, we can still compute balloons shapes at small volumes 
with £>24 symmetry. However, when V « Vj , it may not be reasonable to continue to assume 
that D 24 symmetry is preserved, because of the bulk of excess balloon material hanging at the base 
of the balloon will prevent this from happening. Furthermore, there may be many other factors that 
influence symmetry breaking, load tape slackness (length variation), gore size variation, folding of 
the balloon stack, and other variations introduced while the balloon is in the launch spool. 

2.2 Mathematical Model 

In the numerical calculations presented here on strained lobed balloon shapes, we will use the 
expression for the total energy Ejotal as given by (1.54) and Co = 0 (no aerodynamic loading). 
Furthermore, we will make the assumption that v m [ = v m r = v m = v mr . This will significantly 
reduce the number of degrees of freedom, but we can still obtain reasonable estimates for the 
maximum film stresses 


2.3 Numerical results and discussion 

Figures 2.1- 2.2 suggest that the total energy of a balloon shape is essentially independent of the 
number of lobes when 0.5 V / Vj < 1 . For shapes near V / Vj = 0.2, it appears that shapes with 6 
or 12 lobes have a total energy about 3.5% less than the other shapes for that same volume (see 
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E Total 



Figure 2.1: Total energy for lobed shapes at V/Vj = 0.1, 0.2, 0.3, 0.4, 0.5 (Units of energy lbf-in. 

Figure 2.1). Geometrically, the shapes Sf, and S\2 are very similar. In Sf, there are two flat wings 
per lobe, leading to a total of 1 2 lobes. There is a similar pairing of the shapes, 54 with 5g and S 3 
with Sf,. 

Remark These results suggest that it might be relatively easy for a strained partially inflated 
lobed zero-pressure balloon to transition from a state with n-lobes to one with n + k lobes. It would 
be of interest to study the stability of these shapes and to carry out a similar study for partially 
inflated pumpkin balloons. This could shed some light on problems associated with the ascent and 
deployment of pumpkin balloons. 
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Figure 2.3: Lobed balloon shapes at V/Vj = 0.1, 0.2, 0.3, . . . ,0.9,0.95. (a) Film strain energy; (b) Tape 
strain energy; (c) Hydrostatic pressure potential; 
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(a) o 1 



(b)o 2 



Figure 2.4: Lobed balloon shapes at V/Vj = 0.1, 0.2, 0.3,..., 0.9, 0.95. (a) Max principal stress resultant 
- a i (lbf/in); (b) Max principal stress resultant - q (lbf/in). 
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Figure 2.5: 24-gore balloon: 24 lobes/1 gore per lobe 
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Figure 2.6: 24 


-gore bailoon: 12 iobes/2 


gores per lobe 
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Figure 2.7: 24-gore balloon: 8 lobes/3 gores per lobe 
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Figure 2.8: 24-gore balloon: 6 lobes/4 gores per lobe 
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Figure 2.9: 24-gore balloon: 4 lobes/6 gores per lobe 
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figure 2.10: 24-gore balloon; 3 Jobes/8 


gores per lobe 
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Appendix A 

Aerodynamic drag 


In general, the drag on a submersed body exerted by a fluid that is in motion through the fluid is 


F D — 2 PairC[)A[) | Vg V air |(vg V a(r ), 


(A.l) 


where Co is the drag coefficient, p ai > is the density of the air, Ad is the area of the projection of 
the body on the plane perpendicular to the direction of motion, vg is the velocity of the balloon, 
y a ir is the velocity of the atmosphere (see [7, p. 11-13] and [10, Sec. 4.3]). In our applications, the 
velocity of the balloon will be in the ±k directions. 

Aerodynamic drag can be introduced into our balloon model by adapting Eq. (A.l) appropri- 
ately. Consider an arbitrary facet in the triangulation of a balloon shape that is not part of a fold. 
Let vo, vi , V 2 denote the vertices, and ?o = V] — vo, e\ = V 2 — vj, ei = vo — V 2 . We assume that 
v 0 , vi , V 2 are ordered in such a way that the unit normal n = ?o x e \ /|«o x I points outward. For 
a balloon descending with a velocity vg = — vgk, the drag on typical triangle T is: 



2p ai>CoVg dA • k, 
0 


0 < arccos(— k • n) < 7t/2, 
else, 


(A.2) 


where dA = j|<?o x e\ |n. Adding up the drag contributions of each facet that contacts the atmo- 
sphere yields the total potential energy of the drag on the balloon, i.e., 


E D (S) = ~^F D (T e >), (A.3) 

t' 

where the summation in l 1 is over all facets whose normal satisfies 0 < arccos(— k-n) < Jt/2. 
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Appendix B 

Buoyancy of lifting gas 


In this appendix, we summarize the background related to the calculation of the buoyancy b of the 
lifting gas that is used for the hydrostatic pressure p. While there are numerous parameters and 
combinations of those parameters that lead to the same b, we will provide some of the basic results 
here. We follow the conventions in [7] and [9]. 

Archimedes Principle states that a body submerged in a fluid exerts a force equal to the weight 
of the displaced fluid. In a balloon, the buoyant force or lift is 

B = W a — W g , (B.l) 

where W a is the weight of the displaced air and W g is the weight of the lifting gas. Furthermore, 
we have W a = gp a Vg and W g = g p a V g where g is acceleration due to gravity, p a is the density of 
the air, p g is the density of the lifting gas, and V g is the volume of the gas. The volume of the gas 
is nearly the same as the volume of the displaced air. In particular, we have 


B = gVg(p a -pg) 


(B.2) 


or equivalently 

From the equation of state, we have 

_ PM 
P _ ~RT' 


(B.3) 

(B.4) 


where P is pressure, M is molecular weight, T is absolute temperature, and R is the universal gas 
constant. Using a subscript of V or 'g' for quantities related to air or gas, respectively, we find 


B = gV g p a (\ 


PJaMg\ 

PaT g Mj ‘ 


(B.5) 
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If we assume that 


Eq. (B.5) simplifies to 


T g T a 




In ballooning, the terms superpressure and supertemperature are often introduced (see [7, Section 
II]). Superpressure is defined to be 

n = p g -p a 


and supertemperature is defined to be 


e = n-T a . 


In this case, Eq. (B.5) can be expressed in the form. 


B = gV g Pa 1- 


iPa±ny^Mg\ 

Pa(T a + Q)M a J- 


To determine the buoyancy b of a unit volume of gas in air, we can divide B by Vg and obtain, 


b = gPa(\~^- 


or the more general expression. 


b = gpa ( 1 - 


(p a +n)r a M g \ 
p a (T a + e)M a J' 


(B.10) 


From Eq. (B.10), we see that b = b(p a , r a ,n,0,Af g ,M fl ), 
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Appendix C 

U.S. Standard Atmosphere 


The U. S. standard atmosphere layers are presented in [9, pg. 15]. Data is available on eight 
layers. For each layer, a base layer altitude and top layer altitude is given with corresponding 
values for temperature, pressure, and density. Properties and assumptions are listed in [9, pg. 7]. 
In the following, H is the geopotential height and Z is the geometric height. Within each layer, the 
temperature is assumed to vary linearly from the base of the layer to the top of the layer, i.e., 

T a = T B , i + L' i (H-H B>h Hi<H <H i+u (C.l) 


where Tg,i is the temperature at the base of the ith layer, H, is the base of layer i (the top of layer 
i — 1), and L\ = dT /dH in layer i. From [9, Eq. (12)], we have dP a = —gop a dH 


dH = - 


R 7a 

goM a P a 


dPa- 


From Eq. (C. 1 ) and Eq. (C.2), we have 


dH R dP a 

which can be integrated to yield the height as a function of P a , i.e.. 


(C.2) 


(C.3) 


H=l 


tj , RTb , Pb 

HbI Tr ,n 7T’ 

goM a Pa 


H b + 


h 

L' 


-) 

Pa) 


(L'R)/(g 0 M a ) 


The density of air p a is a function of geopotential height where 


L' = 0, 
L! ^ 0. 


(C.4) 


P a(H) 


PajHWa 
RT a (H) ’ 


(C.5) 
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Table C.l: Physical constants for standard atmosphere 


Quantity 

Sea level air density 
Sea level temperature 
Sea level pressure 
Sea level molecular weight 
Sea level molecular helium 
Sea level molecular hydrogen 
Sea level molecular ammonia 
Universal gas constant 


Value 

1.2250 kg/rtf 
15° C = 288. 15° K 
9.80665 m/sec 2 
28.9644 kg/kg-mol 
28.9644 kg/kg-mol 
2.0159 kg/kg-mol 
4.0026 kg/kg-mol 
8314.32 J/°K(kg-mol) 


and P a ,H, and T a are related by Eqs. (C.l) and (C.4). Related parameters needed to calculate 
atmospheric properties are given in Table C.l. Assuming a standard atmosphere, the buoyancy is 
given by Eq. (B.9) (or Eq. (B.10)), where p a (//) is given by Eq. (C.5), T a is given by Eq. (C.l), 
and H and P a are related by Eq. (C.4). 
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